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Introduction 


For interplanetary missions, highly efficient electric propulsion systems can be 
used to increase the mass delivered to the destination and/or reduce the trip time over 
typical chemical propulsion systems. 1,2 This technology is being demonstrated on the 
Deep Space 1 mission 3 - part of NASA’s New Millennium Program validating 
technologies which can lower the cost and risk and enhance the performance of future 
missions. With the successful demonstration on Deep Space 1 , future missions can 
consider electric propulsion as a viable propulsion option. 

Electric propulsion systems, while highly efficient, produce only a small amount 
of thrust. As a result, the engines operate during a significant fraction of the trajectory. 
This characteristic makes it much more difficult to find optimal trajectories. The 
methods for optimizing low-thrust trajectories are typically categorized as either indirect 
or direct. Indirect methods are based on calculus of variations, resulting in a two-point 
boundary value problem that is solved by satisfying terminal constraints and targeting 
conditions. 4 These methods are subject to extreme sensitivity to the initial guess of the 
variables - some of which are not physically intuitive. Adding a gravity assist to the 
trajectory compounds the sensitivity. Direct methods parameterize the problem and use 
nonlinear programming techniques to optimize an objective function by adjusting a set of 
variables. A variety of methods of this type have been examined with varying 
results. 5,6,7,8 These methods are subject to the limitations of the nonlinear programming 
techniques. 

In this paper we present a direct method intended to be used primarily for 
preliminary design of low-thrust interplanetary trajectories, including those with multiple 
gravity assists. Preliminary design implies a willingness to accept limited accuracy to 
achieve an efficient algorithm that executes quickly. 
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Optimization 


This structure results in a constrained, nonlinear optimization problem which we 
solve using the nonlinear programming software SNOPT 9 The potential set of 
independent variables includes the state (position, velocity, and mass) of the spacecraft at 
each control node and the corresponding epoch. If a control node is associated with a 
solar system body, the position of the spacecraft is the same as the body and therefore is 
not independent. Two independent variables are available for an intermediate flyby, and 
there are variables representing the impulsive AVs on the segments. Depending on the 
optimization objective function and engine model, the solar array reference power and 
engine specific impulse can be independent variables. We normally try to maximize final 
spacecraft mass or net mass (final spacecraft mass - propulsion system mass), but other 
objective functions are possible. 


The primary constraint on the optimization is that the position, velocity, and mass 
of the spacecraft must be continuous at the match points. The magnitude of the impulsive 
AVs may be constrained, as previously described, and other constraints can be placed on 
the trajectory such as total flight time and total propellant mass. In addition, upper and 
lower bounds can be placed on any of the independent variables. 


Results 


We have used the method described in this paper to optimize several different 
types of trajectories. We compare the results to those from SEPTOP, a low-thrust 
trajectory optimization program using an indirect method. (SEPTOP was used in e 
and is briefly described there.) The only major difference between SEPTOP and its well- 
known predecessor VARITOP'^is the way in which the engines are modeled. Both 
programs are the result of a long evolution of low-thrust trajectory optimization software 
and have been used extensively to design a variety of missions. 


We present results from the following three missions in this paper, a flyby of 
Vesta with a Mars gravity assist, a rendezvous with Tempel 1, and a flyby of Pluto with 
two gravity assists at Venus and one at Jupiter. The engine models are based on the 
NSTAR 30-cm ion thruster, 10 a version of which is being flown on Deep Space 1 . e 

initial guesses for the thrust direction and magnitude are crude but simple and have 
worked well. The direction varies linearly between nodes with the direction at the nodes 
being perpendicular to the radius vector of the nodes and in the ecliptic plane The 
maximum AV that could be accomplished on the first segment at a distance ot AU from 
the Sun is used as the starting guess for all segments. 


For the Earth-Mars-Vesta mission, we started out by fixing the launch V„ and 
control node epochs and then subsequently released those variables in a series of runs. 
The final masses from our method and from SEPTOP are shown in Table 1. As can be 
seen from the table, the final masses agree very closely - well within the accuracy of 
either method. Our method converged readily; however, SEPTOP did have some trouble 



Theoretically, the indirect method used by SEPTOP can incorporate an unlimited 
number of intermediate body flybys; however, because of practical limitations arising 
from the sensitivity issues, SEPTOP haS been programmed to accommodate at most two 
intermediate body flybys. Hence, SEPTOP cannot optimize the Earth-Venus-Venus- 
Jupiter-Pluto trajectory in its entirety. To examine such a trajectory, SEPTOP is used to 
optimize the trajectory to Jupiter, and the Jupiter-Pluto leg is determined by C 3 matching. 

Our method can handle any reasonable number of flybys. In fact, using match 
points as a part of the trajectory structure is intended to reduce the sensitivity to adding 
intermediate flybys. To compare to SEPTOP, we fixed the Jupiter flyby and Pluto arrival 
dates, obtaining a final mass of 879.9 kg. Using the procedure described above with 
SEPTOP resulted in a final mass of 880.4 kg. 


Conclusion 

We have developed and tested a direct method for preliminary design of low- 
thrust interplanetary trajectories. This method has been compared to a program using an 
indirect method, and the results agree very closely. The new method has shown less 
convergence sensitivity and the ability to handle more intermediate flybys than the 
indirect method. 
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